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Abstract 

In this article we consider the approximate description of doubly-dispersive channels by its sym- 
bol. We focus on channel operators with compactly supported spreading, which are widely used to 
represent fast fading multipath communication channels. The concept of approximate eigenstructure 
is introduced, which measures the accuracy E p of the approximation of the channel operation as a 
pure multiplication in a given £ p -norm. Two variants of such an approximate Weyl symbol calculus 
are studied, which have important applications in several models for time-varying mobile channels. 
Typically, such channels have random spreading functions (inverse Weyl transform) defined on a 
common support U of finite non-zero size such that approximate eigenstructure has to be measured 
with respect to certain norms of the spreading process. We derive several explicit relations to the 
size \U\ of the support. We show that the characterization of the ratio of E p to some £ g -norm of the 
spreading function is related to weighted norms of ambiguity and Wigner functions. We present the 
connection to localization operators and give new bounds on the ability of localization of ambiguity 
functions and Wigner functions in U. Our analysis generalizes and improves recent results for the 
case p = 2 and q = 1. 

Index Terms 

Doubly-dispersive channels, time-varying channels, Weyl calculus, Wigner function, ambiguity 
function 

I. Introduction 

Optimal signaling through linear time-varying (LTV) channels is a challenging task for 
future communication systems. For a particular realization of the time-varying channel 
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operator the transmitter and receiver design, which avoids interference is related to "eigen- 
signaling". Eigen-signaling simplifies much of the information theoretic treatment of commu- 
nication in dispersive channels. However, it is well known that for an ensemble of channels, 
which are dispersive in time and frequency such a joint diagonalization can not be achieved 
because the eigen-decompositions can differ from one to another channel realization. Several 
approaches like for example the "basis expansion model" (BEM) [1] and the canonical 
channel representation [2] are proposed to describe eigen-signaling in some approximate 
sense. Then a necessary prerequisite is the characterization of remaining approximation errors. 

A typical scenario commonly encountered in wireless communication, is signaling through 
a random time-varying and frequency selective (doubly-dispersive) channel, which in general 
is represented by a pseudo-differential operator 7i. The abstract random channel operating 
on an input signal s : M. — > C can be expressed (at least in the weak sense) in the form of a 
random kernel, symbol or spreading function. The signal r : K — > C at the time instant t at 
the output of the time-varying channel is then: 

r(t) = (Hs)(t) 

It is a common assumption that knowledge of 7i at the receiver can be obtained up to certain 
accuracy by channel estimation, which will allow for coherent detection. However, channel 
knowledge at the transmitter simplifies equalization and detection complexity at the receiver 
and can increase the link performance. It can be used to perform a diagonalizing operation (i.e. 
eigen-signaling) and allocation of resources in this domain (e.g. power allocation). We shall 
call the first part of this description from now on as the eigenstructure of 7i. Signaling through 
classes of channels having common eigenstructure could be, in principle, interference-free 
and would allow for simple information recovering algorithms based on the received signal 
r(t). However, for 7i being random, random eigenstructure has to be expected in general such 
that the design of the transmitter and the receiver has to be performed jointly for ensembles of 
channels having different eigenstructures. Nevertheless, interference then can not be avoided 
in the communication chain. For such interference scenarios it is important to have bounds 
on the distortion of a particular selected signaling scheme. Refer for example to [3] for a 
recent application in information theory. 

Initial results in this field can be found in the literature on pseudo-differential operators 
[4], [5] where the overall operator was split up into a main part to be studied and a "small" 
operator to be controlled. More recent results with direct application to time-varying channels 
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were obtained by Kozek [6], [7] and Matz [8] which resemble the notion of underspread 
channels. They investigated the approximate symbol calculus of pseudo-differential operators 
in this context and derived bounds for the £ 2 -norm of the distortion which follow from 
the approximate product rule in terms of Weyl symbols. We will present more details on 
this approach in Section III-C. Controlling this approximation intimately scales with the 
"size" of the spreading of the contributing channel operators. For operators with compactly 
supported spreading such a scale is \U\ - the size of the spreading support U . Interestingly this 
approximation behavior breaks down in their framework at a certain critical size. Channels 
below this critical size are called according to their terminology underspread and otherwise 
overspread. However, we found that previous bounds can be improved and generalized in 
several directions by considering the problem of approximate eigenstructure from another 
perspective, namely investigating directly the £ p -norm E p of the error Tis — Ar for well 
known choices of A. We shall focus on the case where A is the symbol of the operator 7i 
and on the important case where A is the orthogonal distortion which can be understood as 
the £ 2 -minimizer. We believe that extensions to p ^ 2 are important when further statistical 
properties of the spreading process of the random channel operator are at hand 1 . Our approach 
will also show the connection to well known fidelity and localization criteria related to pulse 
design [6], [9], [10]. In particular, the latter is also related to the notion of localization 
operators [11]. The underspread property of doubly-dispersive channels occurs also in the 
context of channel measurement and identification [12]. In addition refer to the following 
recent articles [13], [14] for rigorous treatments of channel identification based on Gabor 
(Weyl-Heisenberg) frame theory. The authors connect the critical time-frequency sampling 
density immanent in this theory to the stability of the channel measurement. A relation 
between these different notions of underspreadness has to be expected but is beyond the 
scope of this paper. 

The paper is organized as follows: In Section II we shall give an introduction into the 
basics from time-frequency analysis including the Weyl correspondence and the spreading 
representation of doubly-dispersive channels. In Section III of the paper we shall consider 
the problem of approximate eigenstructure for operators with spreading functions, which 
are supported on a common set U in the time-frequency plane having non-zero and finite 
Lebesgue measure \U\. We present the approach for E 2 followed by a summary of the main 

'We provide further motivation and arguments in Remark 2 at the end of the paper. 
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results of our analysis on E, p . The detailed analysis for E p will be presented in Section IV. 
Finally, Section V contains a numerical verification of our results. 

A. Notation and Some Definitions 

We present certain notation and definitions that shall be used through the paper. For 1 < 
p < oo and functions / : R n -> C the functionals ||/|| p := (/ \f{t)\ p dt) lv are then usual 
notion of p-norms (dt is the Lebesgue measure on W 1 ). Furthermore for p = oo is ||/||oo : = 
ess sup If \\f\\p is finite / is said to be in £ p (R n ). The inner product (-, •) on the Hilbert 

space £ 2 (^ n ) is given as (x,y) := L n x(t)y(t)dt where x(t) denotes complex conjugate of 
x(i). A particular dense subset of Cp(R n ) is the class of Schwartz functions 5(IR n ) (infinite 
differentiable rapidly decreasing functions). The notation p' denotes always the dual index 
of p, i.e. 1/p + 1/p' — 1 with p' = oo if p — 1 (and the reverse). 

II. Time-Frequency Analysis 

A. Phase Space Displacements and Ambiguity Functions 

Several physical properties of time-varying channels (like delay and Doppler spread) are 
in general related to a time-frequency view on operators 7i. Time-frequency representations 
itself are important tools in signal analysis, physics and many other scientific areas. Among 
them are the Woodward cross ambiguity function [15] and the Wigner distribution. Ambi- 
guity functions can be understood as inner product representations of time-frequency shift 
operators. More generally, a displacement (or shift) operator for functions / : W 1 — > C can 
be defined as: 

(5^)(x):=e i2 ^V(x-/ii) (1) 

where fi = (/ii,/i 2 ) G M 2n and fii, [12 G W 1 . In general M 2n is called phase space. Later on 
we shall focus on n = 1, where we have that the functions / are signals in time and \i is 
a displacement in time and frequency. Then the phase space is also called time-frequency 
plane and the operators S M are time-frequency shift operators. There is an ambiguity as 
to which displacement should be performed first where (1) corresponds to the separation 
S/j, = S(p ^Sfaoy However, it is well known that a generalized view can be achieved by 
considering so-called a-generalized displacements: 

c(a) — , Cf, . tj _ _-i27r(l/2-a)C(M>/0 <? n\ 

s ti •- a (0,/a(i+a)) i5 ^i ) 0) a (0, w (i-a)) ~ e M V-> 
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where = \i\ ■ v-i (inner product on R n ) and then set = S^ 2 \ Usually a is 

called polarization. The operators in (2) act isometrically on all £ p (IR n ), hence are unitary 
on £ 2 (K n )- Furthermore, they establish 2 unitary representations (Schrodinger representation) 
of the Weyl-Heisenberg group on £ 2 CK ra ) (see for example [5]). In physics it is common to 
choose the most symmetric case a = and the operators are usually called Weyl operators 
or Glauber displacement operators. If we define the symplectic form as 77 (/x, := v) — 
11), we have the following well known Weyl commutation relation: 

g{°<-)g{P) — e -^r)Uhv) g(fi) g(pt) (3) 

for arbitrary polarizations a and (3. In this way a generalized (cross) ambiguity function can 
be defined as: 

A£>0*) d =l f ( 9 , S W 7 ) = f g{x + {\- a)^) 7 (x - (~ + a)^)e l2 ^dx (4) 

The function Ag^ is also known as the Short-time Fourier transform (sometimes also 
windowed Fourier transform or Fourier-Wigner transform) of g with respect to a window 
7. This function is continuous for g E S(R n ) and 7 E iS>'(IR n ) (the dual of i.e. the 

tempered distributions). Well known relations of these functions, which follow directly from 
definition (4) are: 

|A#(A*)I = \{g,Sf 7 )\ < IMI2INI2 = ||A£>|| a (5) 

where the right hand side (rhs) is sometimes also called the radar uncertainty principle. For 
particular weight functions m : R 2n — * R + the weighted p-norms ||Ap" m\\ p are also called 
the modulation norms ||7||m^ p °f 7 with respect to Schwartz function g E S(R. n ) (M^ p is 
then corresponding modulation space [16]). Let the symplectic Fourier transform T S F of a 
function F : IR 2 ™ — * C be defined as: 

{F.F){n)= [ e- l2 ^ v ^F{y)dv (6) 

JR 2n 

The symplectic Fourier transform of the (cross) ambiguity function JF S A^ is called the 
(cross) Wigner distribution of g and 7 in polarization a. 

2 up to unitary equivalence 
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B. Weyl Correspondence and Spreading Representation 

The operational meaning of pseudo-differential operators can be stated with a (distribu- 
tional) kernel, coordinate-based in the form of infinite matrices or in some algebraic manner 
(see for example [17, Chapter 14]). The kernel based description is usually written in a form 
like: 

{Hi){t)= [ h^t'^dt' (7) 

with a kernel h : M 2n — > C (for two Schwartz functions 7, p E 5(IR n ) the kernel h 
exists even as a tempered distribution, i.e. Schwartz kernel theorem states h E S'(M. 2n ) 
with (g,'H'y) = (h,g ® 7), see for example [17, Thm. 14.3.4]). However, the abstract 
description of 7i as superpositions of time-frequency shifts is important and quite close to the 
physical modeling of time-varying channels. We will adopt this time-frequency framework 
to describe the channel operators. Let us denote with the set of compact operators, i.e. 
for X E holds X = J2 k Sk{xk, -)yk with singular values {sk} and two orfhonormal bases 
(singular functions) {xk} and {yk}- For p < 00 the pth Schatten class is the set of operators 
T p := {X \ \\X\\p := Tr((X*Xyl 2 ) = J2 k \ s k\ P < °°} wh ere Tr(-) is the usual meaning 
of the trace (e.g. evaluated in a particular basis). Then T p for 1 < p < 00 are Banach 
spaces and T x C T p C (see for example [18]). The sets 7[ and T 2 are called trace class 
and Hilbert-Schmidt operators. Hilbert-Schmidt operators form itself a Hilbert space with 
inner product (Y,X)t 2 := Tr(Y*X). For X E T\ it holds by properties of the trace that 
\(Y,X)r 2 \ < ||X||i||y||, where ||-|| denotes the operator norm. Hence for Y = given 
by (2) one can define analogously to the ordinary Fourier transform [19], [20] a mapping 
7! -> £ 2 (M 2n ) via: 

vP(p)^(S<?\X)t 2 (8) 

In essence, the kernel h of the channel operator 7i is given as the (inverse) Fourier trans- 
form in the p 2 variable (see for example [17, Chapter 14]). Note that £^(0) = Tr(X) 
and < ||X||i (because HS^H = 1). The function is sometimes called the 

"non-commutative" Fourier transform [21], characteristic function, inverse Weyl transform 
[22] or a— generalized spreading function of X [6], [23]. From (2) it follows that — 

e i27r(l/2-a)C(M,A05] (1/2) . 

Lemma 1 (Spreading Representation) Let X E T 2 . Then there it holds: 

X= f (S^\X) T2 S^dp (9) 

JR 2n 
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where the integral is meant in the weak sense 3 . 

The extension to the Hilbert-Schmidt operators T 2 is due to continuity of the mapping in (8) 
and density of T x in T 2 . A complete proof of this lemma can be found in many books on Weyl 
calculus (for example matched to our notation in [21, Chapter V] and [24]). Furthermore, 
the following important shift-property: 

can be verified easily using (2) and (3). The composition of the symplectic Fourier transform 
T s as defined in (6) with the mapping in (8) establishes the so called Weyl correspondence [22] 
in a particular polarization a (for this generalized approach in signal processing see also [23]). 
The function = JF S £^ is called (generalized) Weyl symbol of X. The original Weyl 
symbol is The cases a = \ and a — — \ are also known as Kohn-Nirenberg symbol 
(or Zadeh's time-varying transfer function) and Bello's frequency-dependent modulation 
function [25]. The Parseval identities are: 

(X,Y) T2 = = <4 a U? } > (ID 

for X, Y e T 2 . For a rank-one operator X = (7, -)g it follows that = A^y such that 
(11) reads in this case as: (g,Y^y) = (A^ , £y '). 

III. Problem Statement and Main Results 

In this section we will establish a concept, which we have called the "approximate eigen- 
structure". The latter are sets of signals and coefficients which fulfill a particular property of 
singular values and functions up to certain approximation error E p measured in p-norm. Part 
III-A motivates this concept for a single channel operator. In part III-B of this section we 
will then extent this framework to a time-frequency formulation for "random" time-varying 
channels with a common support of the spreading functions. We consider on how approximate 
eigenstructure behavior scales with the respect to the particular spreading functions, which 
is the main problem of this paper. Recent results in this direction are for p = 2 and based 
on estimates on the approximate product rule of Weyl symbols. We will give a general 
formulation of this approach and an overview over the known results for E 2 in part III-C 
of this section. After that we present in III-D a new (direct) approach for upperbounding 

3 For HXj^' ||i < 00 (9) is a Bochner integral. Weak interpretation of (9) as (g,Xj) extents the meaning of this integral 
to tempered distributions [5, Chapter 2] or [17, Chapter 14.3]. 
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E p yielding for our setup also improved and more general estimates for p — 2. This part 
contains a summary of the main results, where the more detailed analysis is in Section IV. 

A. The Approximate Eigenstructure 

It is a common approach to describe a given channel operator 7i on a superposition of 
signals where 7i act rather simple. As already mentioned, compact operators on a Hilbert 
space can be formally represented as 7i = Y^k=i s k{ x k, -)Vk with the singular values {s*,} 
and singular functions {xk} and {yk}- Transmitting an information bearing complex data 
symbol c for example in the form of the signal s = c ■ Xk through H, we known that with 
proper channel measurement (obtaining Sk) the information can be coherently "recovered" 
from the estimate (yk, 7~Cs) = s& • c. The crucial point here is that the transmitting device has 
to know and implement {xk} before. However, in practical implementation is required 
to be fixed and structured to some sense (for example in the form of filterbanks). But in 
general, also the singular functions depend explicitly on the operator 7i, i.e. they vary from 
one realization to another. They can be very unstructured and it is difficult to relate properties 
of 7i in such representations to physical measurable quantities. 

Hence, instead of requiring Tixk = SkVk we would like to have that Tixk — Sf-Vk is "small" 
in some sense. Usually, approximation in the £ 2 -norm seems to be of most interest in the 
signal design. However, there are certain problems as peak power and stability issues where 
stronger results are required. Furthermore intuitively we are aware that the approximation of 
the singular behavior of {s*., yk, x^} has to be "uniform" in more than one particular norm. 
In this paper we consider the C p norms for the approximation, thus we have the following 
formulation for the Hilbert space £ 2 (^ n ): 

Definition 2 (Approximate Eigenstructure) Let e be a given positive number. Consider 
A G C and two functions gjG /^(K™) with \\g\\2 = IMh = 1. If 

E p := \\Hj-Xg\\p<e (12) 

we call {A, g,j} a C p -approximate eigenstructure ofH with bound e. 

The set of A's for which exists g\ such that {A, g\, g\} is a £ 2 - a Pproximate eigenstructure for 
a common fixed e is also called the e-pseudospectrum 4 of 7t. More generally, we will allow 

4 Thanks to T. Strohmer for informing me about this relation. 
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also g ^ 7 such that the term "approximate singular" functions is suited for our approach 
as well. Obviously for the "true eigenstructure" {sk,yk,Xk} as defined above we have that 
e = for each p and k . On the other hand, for given g and 7 the minimum of the left hand 
side (lhs) of (12) is achieved for p = 2 at A = (gjHry) such that E p for {{g, Hr/), g, 7} 
describes the amount of orthogonal distortion caused by TL measured in the p-norm. 

B. The Problem Statement for Channels with Compactly Supported Spreading 

It is of general importance to what degree the Weyl symbol or a smoothed version of 
it approaches the eigen-value (or more generally singular value) characteristics of a given 
channel operator 7i. Inspired from the ideas in [7] we will consider now the following ques- 
tion: What is the error E p {p) if we approximate the action of 7i on S^j as a multiplication 
of S^g by A(/x)? Hence, instead of the "true" eigenstructure consisting of the singular values 
and functions of 7i we shall consider a more structured family {A(/x), S^g, S^}. The latter 
will intuitively probe the operator 7i locally in a phase space (time-frequency) meaning if 
g and 7 are in some sense time-frequency localized around the origin. The validity of this 
approximate picture, in which the function A : M 2n — > C now serves as a multiplicative 
channel is essentially described by E p (p). 

For example, in wireless communication S^g and S^j could be well-localized prototype 
filters at time-frequency slot p of the receive and transmit filterbanks of a particular com- 
munication device and \(p) is an effective channel coefficient to be equalized. However, 
with this application in mind, we are typically confronted with random channel operators 7i 
characterized by random spreading functions £^ having a common (Lebesgue measurable) 
support U of non-zero and finite measure \U\, i.e. < \U\ < 00. The assumption of a known 
support seems to be the minimal apriori channel knowledge that enters practically the system 
design (e.g. of a communication device). For example, a typical doubly-dispersive channel 
model (n = 1) for this application is that spreading occurs in U = [0,r d ] x [—B D ,B D ] 
where and -Bd are the maximum delay spread and Doppler frequency. It is then desirable 
to have common prototype filters for all these channel realizations. It is clear that in this 
direction Definition 2 is not yet adequate enough. We have to measure the approximation 
error with respect to a certain scale of the particular random spreading functions. In this 
paper we measure the approximate eigenstructure with respect to its £ g -norm. We believe 
that this approach is important to have reasonable estimates for the various statistical fading 
and scattering environments. An example of such an application is given in Remark 2 in 
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Section V. 

We consider only bounded spreading functions such that the operator 7i is of Hilbert 
Schmidt type, i.e. 7i G T 2 . To this end, let us call this set of channel operators as OP(U), 
i.e. 

OP(U) := {H |supp(sg } ) C [/and sup < 00} (13) 

As already discussed for example in [14] the operator class OP(£7) does not include limiting 
cases of doubly-dispersive channels like the time-invariant channel or the identity. Gener- 
alizations, for example in the sense of tempered distributions, are beyond the scope of this 
paper. We aim at an extension of Definition 2 for the approximate eigenstructure which is 
meaningful and suited for this class of channels. We will formulate this as our main problem 
of this paper: 

Problem: Consider two functions g, 7 : MJ 1 — > C with \\gW2 = W1W2 = 1- Let be 1 < q < 00, 
1 < p < 00 and < 5 < 00 such that for all operators 7i G OP(U) it holds: 

E^) : = \\HSfi - \<ji)Sr9\\ p < 5 ■ HE^II, (14) 

where the p-norm is with respect to the argument of the function TiS^j — X^S^g. Then 
{A(/z) , S fj,g , S ^7} is an C p -approximate eigenstructure for all H G OP{U), each of them 
with individual bound e = 5 ■ ||£^|| a . How small can we choose the scale 5 given g, 7, U, 
p and q? What can be said about inf 97 (<5)? 

Note that, independently of the polarization a, the operator S p can be replaced in (14) with 
any /5-polarized shift without change of E p (n). Furthermore, as already stated in the 
definition of E p in (12) ||g|| 2 = H7H2 = 1, throughout the rest of the paper. Summarizing: 
How much could {S^g, S^} serve as common approximations (measured in p-norm) to 
the singular functions of the operator class OP(£7) for fixed U ? 

C. Results Based on the Approximate Product Rule 

In previous work [6], [26], [8], [27] results were provided for g = 7 and (apart of [28]) 
A = for the case p = 2. These are obtained if one considers the problem from view of 
symbolic calculus and can be summarized in the following lemma: 

Lemma 3 Let 7 = g and X = L^. It holds: 

em < (iL^M - |££V)I 2 I + lis^lli + 2|x£V)l • IIEgta Ik) 1 (15) 
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where Vt = |A 7 " — 1|. 

Note that the lemma not yet necessarily requires 7i G OP(£7). The proof is provisioned in 
Appendix A. This bound is motivated by the work of W. Kozek [6]. However it has been 
formulated in a more general context. The first term of the bound in (15) contains the Weyl 
symbol L^y of the composition XY of two operators X and Y (7i* and 7i in this case), 
which is the twisted multiplication [29] of the symbols of the operators X and Y. On the 
level of spreading functions 5 , £^ y is given by the so called twisted convolution \ $ of 
and e£° [24], [30]: 

SS(P) = / SS? J W4% - tie-^^df, (JW i, , S W) (p ) (16) 

with 0(p, p) = (a + |)£(p, p) + (a — |)C(p> p) — 2«C(p, p). For the polarization a = 
it follows (f)(fi,p) = p(p,p)/2 and conventional convolution is simply t| . Expanding 
exp(— z27T0(/i, p)) in p as a Taylor series reveals that twisted convolutions are weighted sums 
of \ o-convolutions [5] related to moments of and £y ■* . Hausdorff- Young inequality 
with sharp constants = pp jpp' (and c\ = = 1) gives for 1 < p < 2 estimates on the 
following "approximate product rule" of Weyl symbols: \\L^ Y — L^Ly^Wp' < 2c£ n ||-F||p 
where F{p) = j R2n |£^(p)£y ^(p — p) sin(27T0(p, p))\dp. In particular, for p = 1 we get: 

I^S*- l^5?l 2 l < 2 ll F lli a - e - < 17) 

Let us assume now that Ti E OP(U). With xu we shall denote the characteristic function of 
U (its indicator function). Kozek [6, Thm. 5.6] has considered the case a = obtaining the 
following result: 

Theorem 4 (W. Kozek [6]) Let U = [-t ,t ] x [-z/ ,^o] and a = 0. If \U\ = At u < 1 
then 

E 2 (p)< ( 2sm (^)||E«?||? + e J (\\^ n \\ 1 + 2\\^\\j)y (18) 
where e 7 = ||(A$ - l)xi/||oo- 

The proof can be found in [6] or independently from Lemma 3 with e 7 = || f^Xfy II oo and 
H-k^lloo < ll^w Ik- Further utilizing the fact that HE^l^Hi < ||S^||f, Equation (18) can 

5 Symbols (spreading functions) in £2(R 2n ) with twisted multiplication (convolution) are *-isomorph to the algebra of 
Hilbert-Schmidt operators. 
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be written as: 




CO 



) 



2 



(19) 



which gives an initial answer to the problem formulated in section III-B. Theorem 4 was 
further extended in [8, Thm. 2.22] by G. Matz (see also [27]) to a formulation in terms 
of weighted 1-moments of (not necessarily compactly supported) spreading functions. His 
approach includes also different polarizations a. For a spreading function as in Theorem 4 
and a = the results agree with (19). Equation (19) could be interpreted in such a way that 
only the second term can be controlled by 7 (e.g. pulse shaping) where the first term of the 
rhs of (19) is only related to the overall spread \U\. However we shall show in the next section 
that the first term can be eliminated from the bound and the second (shape-independent) term 
can further tightened. 

D. Results Based on a Direct Approach 

We have considered the function A as the Weyl symbol exclusively for the exponents p = 2 
and a = 1 in III-C. This approach is in line with prior work of Kozek, Matz and Hlawatsch 
and provides results on the approximate eigenstructure problem established in Section III-B. 
To obtain further results for different values of p, a and A, we shall now restart the analysis 
from a different perspective. In the following we present the main results, through most of the 
analysis will be presented in Section IV. We use a "smoothed" version of the Weyl symbol: 



where B : IR 2n — > C is a bounded function. We consider two important cases: 

Case CI: Let B = Ag$ such that (20) reads as A = * JT S A^ where * denotes 
convolution. This corresponds to the well known smoothing with the cross Wigner function 
and was already considered in [28] (for averages over WSSUS 6 channels ensembles). 
In particular this is exactly the orthogonal distortion: 



as already mentioned in III-A and corresponds to the choice of the E 2 -mmvmize,x. Since A 
depends in this case on g and 7 we consider here how accurately the action of operators 7i 
on the family {Syr} can be described as multiplication operators on the family {S^g}. From 

6 Wide-sense stationary uncorrelated scattering (WSSUS) channel model [25] 



A = r 8 (Y% ■ B) 



(20) 



(21) 
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the rule in (2) the definition of the cross ambiguity function in (4) and the non-commutative 
Fourier transform in (8) it is clear that this choice is also independent of the polarization a. 
Recall that the Weyl symbol of a rank-one operator is the Wigner distribution, such that with 
this approach we again effectively compare twisted with ordinary convolution. 

Case C2: Here we consider B = 1 such that A = is the Weyl symbol. The function 
A is now independent of g and 7. Thus in contrast to CI this case is related to the "pure" 
symbol calculus. Obviously, we have to expect now a dependency on the polarization a. 
Furthermore, this was the approach considered for p = 2 in the previous part 7 of this section. 

The first theorem parallels Theorem 4 and its consequence (19). We shall not yet restrict 
ourselves to the cases CI and C2. Instead we only require that B has to be essentially 
bounded. 

Theorem 5 Let H G OP{U), jjG £oc(K") and B e £oc(M 2rt ). For 2 < p < 00 and 
1 < q < 00 (with the usual meaning for q = 00): 

-^L < C ^.\\(1+ |5| 2 - 2Re{A^B}) Xu \\% (22) 

where C is a constant depending on g, 7 and B. The minimum of this bound over B is 
achieved in the case of p = 2 for CI. 

Proof: The proof follows from the middle term of (39) in Lemma 9 if we set C as: 

C = ess sup |(5[f ) 7)(a;) - B(u)g(x)\ < IMU + ||5|Ub||oo (23) 

x£R n , u&J 

In Lemma 9 the range of p is 1 < p < 00. However, from the discussion in Section IV-B it 
is clear that (22) gives only for p > 2 a reasonable bound. ■ 

Comparison to the bound of Kozek: With |1 — Re{A^}| < |1 — aJ$*| we can transform 
the result of the last theorem for C2 with settings p = 2, q — 1 and g = 7 into a form 
comparable to (18) and (19) which is: 

EM 



<(2||(l-A^)xdloo) 5 (24) 



iv»(«)ii ~ v " v 97 
l^n Hi 



Hence this technique improves the previous bounds. It includes different polarizations a and 
does not require any shape or size constraints on U. Interestingly the offset in (19), which 

7 However, also there the same methodology as in (20) could be applied as well. 
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does not depend on (g,j) and in an initial glance seems to be related to the notion of 
underspreadness, has been disappeared now. 

Discussion of the critical size: The behavior of the bound in (22) on \U\ depends in 
general on the choice of the function B. For example for the case CI, p = q = 2 and with 
(5) it follows that the rhs of (22) is the square root of \U\ — (|A^| 2 , \u) and again with (5) 
we have that: 

y/\U\ -min(|C/|,l) < rhs of (22) < y/\U\. (25) 



This implies that this term is of the same order as y/\U\ for \U\ ^> 1 (see also Lemma 
11 later on). The lhs of the inequality suggests that for \U\ < 1 the scaling behavior might 
alter, i.e. \U\ — 1 is in this sense a critical point between over- and underspread channels as 
introduced in [6]. On the other hand the lhs of the last equation is not zero for < \U\ < 1. 
Indeed from Theorem 14 we have an improved version as follows: 



\U\ - min(|[/|e-^, 1) < rhs of (22) < ^\U\ (26) 

which suggest that at \U\ = e the behavior changes. 

Restriction to the cases CI and C2: If we further restrict ourselves to q > 1 (i.e. q' < oo) 
we can establish the relation to weighted norms of ambiguity functions [10]. For simplicity 
let us consider now the two cases CI (k = 1) and C2 (k = 2). We define therefore the 
functions A k : R 2n -> R for k = 1, 2 as: 

Ai:=|A$| 2 and A 2 := Re{A^} (27) 

We then have the following result: 

Theorem 6 Let H e OP{U) and 9,76 C^iW 1 ). For 2 < p < 00, 1 < q < 00 and \U\ < 1 
it holds: 



iv( Q )ii 
l^n Ik 



< C v ~k {k{\U\ - (A k , xu))) 1/max(q '' p) (28) 



where k = 1 for CI and k = 2 for C2. 

Proof: We now use the bound (40) in Lemma 9 with the uniform estimates Cb p < k 
from Lemma 10. Again, as follows from the discussion in Section IV-B, we consider only 

P > 2. ■ 

The assumption \U\ < 1 is only used to simplify the bound. Improved estimates follow from 
Lemma 9 directly. From the positivity of Ai we observe that the orthogonal distortion (the 
case CI) is always related to weighted 2-norms of the cross ambiguity function (the weight 
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is in this case only xu)- For the case C2 this can be turned into weighted 1-norm if A 2 
is positive on U or fulfill certain cancellation properties. Furthermore, the case C2 depends 
obviously on the polarization a. For particular symmetries of U explicit values can be found 
as shown with the following theorem (for simplicity we consider n — 1): 

Theorem 7 If H e OP(U) with xv{p) = Xu{-^) and 9,1 e ^oo(^) and B e £oo(R 2 ). 
For 2 < p < oo, 1 < g < oo and |Z7| < 1 «Y holds: 

Ep ^ <32 P -^k(k\U\(l-L)) 1/miixW ' p) (29) 



In general, L > \ ma , x (Q*Q) 1 ^ k and Q is an operator with spreading function Xu/W\ in 
polarization 0. Furthermore, L > l(\U\2/k) with l(x) = 2(1 — e~ x ^ 2 ) /x for U being a disc 



and l(x) = 2 • erf[\Jnx/%) 2 /x for U being a square. 

Proof We combine Lemma 9 and Lemma 13 with the uniform estimates < k from 
Lemma 10 such that 

< 32^k (k\u\(i - \ m UQ*Q) 1/k )) 1/maxW ' P) OO) 

where Q is a compact operator with spreading function Xu/\U\ in polarization a. From 
Lemma 15 we know that our assumptions imply that Q is Hermitian for a = 0. In general, 
therefore it holds that L = A max (Q*Q) 1 ^ = A max (Q) 2 ^ where A max (Q) is at least as the 
value of the integral (56) over the first Laguerre function. We abbreviate L = l(\U\) 2 ^ h such 
that for a disc of radius \/\U\/n this integral is l(x) = 2(1 — e~ x ^ 2 )/x and for a square of 
length y/TJ we get l(x) = 2 • erf( yJitx/S) 2 / x. However, Lemma 15 asserts this as an upper 
bound achieved in this case with Gaussians g(x) = 2 1 ^e~ n( - x ' x ' ) which is tight for C2 but not 
for CI. This means, for CI this can be further improved by direct evaluation on Gaussians. 
Indeed, from the proof of Corollary 16 we know that l(\U\ ■ 2/k) > l(\U\) 2 ^ k is achievable. 
For Halloo we get \\VWn = 2^^ = 32 1 / 4 . ■ 
This result can be extended in part to regions U which are canonical equivalent to discs and 
squares centered at the origin (see the discussion at the beginning of Section IV-D). However, 
this holds in principle only for p = 2 because such canonical transformations will change 
the constants in (29). 

IV. General Analysis and Proofs 

With the following lemma we separate the support of the spreading function £^ from 
the quantity E p (p). We shall make use of the non-negative function V : W l x R 2n — > R + , 
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defined as: 

V(x, v) := \(Si a h)(x) - B{v)g{x)\ ■ X u{v) > (31) 

and of the functionals V p {v) := \\V(-, v)\\ p , i.e. the usual p-norms in the first argument. For 
simplifying our analysis we shall restrict ourselves to indicator weights xu (the characteristic 
function of U). However, the same can be repeated with slight abuse of notation using more 
general weights. 

Lemma 8 Let U G OP{U), 1 < p < oo, 1 < q < oo. If V(-, v) G C p (R n ) for all v G U 
then it holds: 

EM/\\Z$\\ q < \\V p \\ q , (32) 
whenever G C q (R 2n ) and V p G C q ,(R 2n ). 

Proof: Firstly, using Weyl's commutation rule (3) and the definition of A in (20) gives: 

em ^ || / d^(u)s^s pl - ms,g\\ P 

JR 2n 

( => \\sjf du^\u)e- l2 ^Si a h ~ \(M)g) Hp (33) 
W || f dv s£ } (v)e-^M (5 («) 7 _ B{y)g) \\ p . 

Note that the p-norm is with respect to the argument of the functions g and S^j. The last 
step follows because 5j t o) acts isometrically on all C p (R n ). Let / : R n x R 2n -> C be the 
function defined as: 

f(x, v) := e^sjj)^)^)^) _ B(u)g(x)} (34) 

From H G OP(U) (bounded spreading functions) and V{-,v) G C p (R n ) for all v G U it 
follows that f(-,v) G C p {R a ). Then (33) reads for 1 < p < oo by Minkowski (triangle) 
inequality 

EM = \\[ duf(-,u)\\ p < [ ^/(^IIpHI^-^IIi^IIS^IUI^L' (35) 
In the last step we used Holder's inequality, such that the claim of this lemma follows. ■ 
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A. The Relation to Ambiguity Functions 

In the next lemma we shall show that ||V^|| g / can be related to ambiguity functions, which 
occur for p = 2. We introduce R : M? n — > M + as the non-negative function: 

R := V 2 2 = (1 + |A$ - 5| 2 - |A^| 2 ) • X t/ > (36) 

and abbreviate R s := ||-R|| S . Using (27) we can write (36) for the cases CI and C2 as 
R = k(l — Ak) ■ Xu for k = 1,2. From the non-negativity of R follows that: 

R 1 = \U\ + ||(A( Q 7 ) - B) Xu \\l - \\A$ X u\\l (37) 

Hence, R\ reflects an interplay between two localization criteria in the phase space. In 
particular, we get for CI and for C2: 

R 1 = k(\U\-(A k , Xu )) (38) 

With the following lemma we shall explicitly provide the relation between the bound ||VJ,||g> 
in Lemma 8 to the quantity R x . 

Lemma 9 For 1 < p < oo and 1 < q' < oo it holds (with the usual meaning for q' = oo): 

ll^|| ? '<||^||5 2 -^ (39) 
Equality is achieved for p = 2 and then the minimum over B of the rhs is achieved for CI. 

Q 1 — V I 

- f - p—q 

For q f < oo let C pq = R<£ p for p < q f and C pq = \U\ else. Then it holds further that: 

\\V P \W< \\V\\T ■C pg -R 1 1 / - xM) (40) 

with equality for q' — p — 2. 

The proof can be found in the Appendix B. The main reason for this lemma, in particular for 
the second part, is that it opens up for case CI the relation to weighted norms of ambiguity 
function (i.e. localization of A k on U). However, for C2 we are also concerned with the 
question of positivity (and cancellation properties) in U. We shall study these relations in 
more detail in Section IV-C. 
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B. Uniform Estimates 

As already mentioned before, for "true" eigenstructure we have E p = for all p, such that 
the notion of approximate eigenstructure should be in some sense uniform in q and p. In the 

p-2 

first step it is therefore necessary to validate uniform bounds for C pq . We observe that || V||oo 
will then restrict the application of Lemma 9 only to p > 2 because || V\\oo will be in general 
small. For example for C2 and g = 7 let \U\ — > in (31). This behavior has to be expected 
because the ambiguity function is a £ 2 -rel a ted construction and from £ 2 boundedness one 
can only with further decay conditions infer £ p -boundedness for p < 2. Consequentially we 

p-2 

shall restrict the following analysis to 2 < p < 00 such that sup||\/||oo p = max(|| V||oo, 1). 
For || V || 00 we can use for example a worst case estimate of the form ||V||cc < ||t||oo + 
||-B||oo ■ \\g\\oo < IItIU + Iblloo which is valid for CI (||-B||oc = II A 97 ||oo < 1 by (5)) and C2. 



Lemma 10 (Uniform Bounds for C pq ) For 2 < p < 00 and 1 < q < 00 it holds the 
uniform estimate C pq < k if q' > p where k = 1 for CI and k = 2 for C2. If q' < p then it 
holds C pq < max(|{7|, 1). 

p — q' 

Proof: It is easily verified that sup If/) - ^ = max(|£/|, 1) where the supremum is over 
all 1 < q' < p and 2 < p < 00. The same can be found also for 1 < p < 00. Similarly we get 

q'-p q' —p 

for the quantity R^ p the uniform estimate sup R^ p = max(^/R^ l , 1) where p < q' < 00 
2 < p < 00. For 1 < p < 00 we would get instead max(i? 00 , 1). From the non-negativity of 
R it follows that: 

Roe = k{l - ess inf A k {v)) (41) 

From (5) it follows that the inequality R^ < 1 is always fulfilled for CI. For the case C2 
this gives instead that i?^ < 4, in general. ■ 
The following lemma provides a simple upper bound on £ , p /||S^' ) ||q which is for p = 2 
uniformly in g and 7. Thus, it will serve as a benchmark. 

Lemma 11 (Uniform Bound for £? p /||S^ || g ) For 1 < p < 00 and 1 < q < 00 it holds: 

EM < \\V\\f ■k 2 / p -\U\ 1 ^ (42) 



\ Z -"H II? 



with k = 1 for CI and k = 2 for C2. 

Proof: We use Ri < R^ ■ \U\ in (39) of Lemma 9 and the uniform estimates R^ < k. 
from Lemma 10. ■ 
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This bound can not be related to ambiguity functions, i.e. will give no insight on possible 
improvements due to localization. 

C. Weighted Norms of Ambiguity Functions and Localization 

In the previous section we have shown that Ri is a relevant term, which controls the 
approximate eigenstructure. In the following analysis we shall further investigate R±. We are 
interested in M gr/ (Ri) which is: 

inf i*i = k\U\ (1 - sup(A fe , C) ] (43) 

™ V 9,7 / 

where C := Xu/W\- Thus, (43) is a particular case of a more general problem, where C 
is some arbitrary weight (non-negative) function C. Thus, let us consider sup gr/ (Ak,C) 
and let us focus first only on A\ = lA^I 2 which is also positive. Since A\ is quadratic 
in 7 we can rewrite (Ai,C) = (7, Lc, g l) where this quadratic form defines (weakly) an 
operator Lc, g - Such operators are also called localization operators [11] and it follows that 
sup 7 (Ai,C) = A max (Lc, i ,). The eigen-values and eigen-functions of Gaussian (g is set to 
be a Gaussian) localization operators on the disc (U is a disc) are known to be Hermite 
functions (more generally this holds if C has elliptical symmetry). Kozek [6], [7] found that 
for elliptical symmetry also the joint optimization results in Hermite functions 8 . For C being 
Gaussian the joint optimum (g and 7) is known explicitly [10]. The last result is based on 
a theorem, formulated in [10], which we will need also in this paper. Let us consider for 
simplicity once again the one-dimensional case (the generalizations for n > 1 are similar), 
i.e. for n — 1 we have: 

Theorem 12 Let \\gW2 = 1 1 T 1 1 2 = 1 and s,r e R Furthermore let C E £ S /(IR 2 ). Then the 
inequality: 

<|A^,C0<(|)V|| S ' (44) 

holds for each s > max{l, 2 }. 

From (2) follows that (44) does not depend on the polarization a. The proof can be found 
in [10] and is based on a result of E. Lieb [31]. Note that apart from the normalization 

8 Kozek considered 3 = 7. However one can show that for elliptical symmetry around the origin the optimum has also 
this property. 
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constraint the bound in Theorem 12 does not depend anymore on g and 7. Hence for any 
given C the optimal bound N T (C) can be found by 



1 



N r (C) := min I ( — ) ||C||, I (45) 

-5 S >max{l,2} \\rsj J 

The equality case in Theorem 12 is given for #,7 and C being Gaussians (see [10] for 
more details). The following lemma states lower and upper bounds on the optimal achievable 
values of the quantities (A k , C). 

Lemma 13 Let be C : IR 2n — > R + a non-negative weight function with \\C\\i = 1. Then it 
holds: 

A max (Q*Q) < sup(A 1 , C) < N 2 (C) (46) 

9,7 

for case CI and equivalently for case C2: 

KUQ*Q) l/2 = max(A 2 , C) < N^C) (47) 

9.7 

where Q is the operator with spreading function C in polarization a. 

Proof: Considering first the case CI (that is k = 1), which is independent of the 
polarization a. The corresponding term (Ai,C) is relevant in the theory of WSSUS pulse 
shaping [9] where C is called the scattering function. In [32] we have already pointed out 
that a lower bound can be obtained from convexity. We have: 

|(^,Q7)| 2 <(^i,C)<iV 2 (C) (48) 

where Q is a compact (follows from normalization) operator with spreading function C. The 
uniform upper bound is according to (45). The optimum of the lower bound is achieved for g 
and 7 being the eigen-functions of Q*Q and QQ* corresponding to the maximal eigen-value 
Kmx(Q*Q), such that for the supremum over g and 7 it follows that: 

KUQ*Q) < snp(A 1 , C) < N 2 (C) (49) 

9,7 

For the case C2 (k = 2) we proceed as follows. For a given 7 we have: 

(A 2 , C) = i «Q 7 , 9) + (9, Qi)) < WQlh (50) 

with equality in the last step for g = Qj/\\Qj\\2- Choosing 7 from the eigen-space of Q*Q 
related to the maximal eigen-value, we get: 

A max (Q*Q) 1/2 = max(A 2 , C) < N^C) (51) 

9,7 
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because (A 2 , C) <(\A 2 \,C) < (V^T, C) < N^C) where again N x is from (45). ■ 
For the particular weight function of interest in this paper, i.e. for C = Xu/\U\ the upper 
bounds can be calculated explicitely. For n = 1 we get the following result: 

Corollary 14 (Norm Bounds for Flat Scattering) Let be C := Xu/\U\. Then it holds that: 

{e-^ \U\ < 2e/r* 

(HW else 

where r* = max{r, 2}. It is not possible to achieve equality. 
The proof is obviously independent of a and available in [10]. 

Remark 1 When using the WSSUS model [25] for doubly-dispersive mobile communication 
channels one typically assumes time-frequency scattering within a shape U = [0, r d ] x 
[—Bd, Bd] such that \U\ = IB^Jd C 1 < e, where Bd denotes maximum Doppler bandwidth 
Bd and Td is maximum delay spread. Then (52) predicts for a L\-normalized scattering 
function C := \U\~ 1 Xu> tnat the best (mean) correlation response (r = 2) in using filter g 
at the receiver and 7 at the transmitter is bounded above by e~ 2BdTd ^ e . 

From the definition of R± in (37) and from (52) of Corollary 14 we know that for \U\ < ke 
we have the estimate: 

k\U\(l - e"^) < inf^) < k\U\(l - \ m UQ*Q) 1/k ) (53) 

9,1 

which are implicit inequalities for \U\. The restriction \U\ < e for the lower bound can be 
removed if the second alternative in (52) of Corollary 14 is further studied. However, for 
simplicity we have considered only the first region which is suited to our application (small 
\U\). In particular, with Ri < Roo\U\ we have also R^ > k(l — e~^). This proves also the 
assertion in [33], i.e. a necessary condition for R^ < 1 is that \U\ < 2eln2 . Furthermore 
for i?oo — > k the size constraint on U vanishes. 

D. Even Spreading Functions and Laguerre Integrals 

Simple estimates for (|A^y| r >C) (and therefore also for X max (Q*Q)) can be found if 
C exhibits certain symmetries upon canonical transformations. Let T : M. 2n M 2n be the 
transformation T(v) = L ■ v + c with a 2n x 2n symplectic matrix 9 L and a phase space 

9 This means that r](L/i, Lfi) = r/(n, fi) for all /i. In particular this means that | det(L)| = 1 such that the measure \U\ 
is invariant under L. 
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translation c £ K . It is well known that |A 97 | = |Ag 7 o T\, where g and 7 are related 
to g and 7 by unitary transforms which depend on T. See for example [5, Chapter 4] for 
a review on metaplectic representation. We have then: (|A 97 | r ,C) = (|Ag 7 | r ,C 0T" 1 ). In 
particular this means, that we can always rotate, translate and (jointly) scale C to simple 
prototype shapes. For example, elliptical (rectangular) shapes can always be transformed 
to discs (squares) centered at the origin. Further symmetries can be exploited as shown 
exemplary in the following lemma (for simplicity we consider only n — 1): 

Lemma 15 Let be Q the operator with spreading function xu- If the shape of U has the 
symmetry Xu{fJ>) = Xu{—^) then for each m > it holds that: 

A m ax(<5*Q) > (tttt / l m {*W))dy\ (54) 



\\U\ju 

where |/i| 2 = /i 2 + /i 2 . and l m is the mth Laguerre function. 

Proof: The calculation of X max (Q*Q) simplifies much for normal operators which 
involves the investigation of Q only, i.e. X m (Q*Q) = |A m (Q)| 2 . For an arbitrary operator Y 
it follows that £y](/i) = s|T (— ^e~ li7ra ^'^ is the spreading function of Y* in polarization 
a. Hence, on the level of spreading functions the normality of Y is equivalent to: 

E^V)S5?V) = T^\-^\-v) ■ e u*««M+«w)) (55) 

which can be verified using the rules for like (2) and (3). The operator Q has by 
definition the real spreading function xu- Hence the desired symmetry is fulfilled for a = 0. 
Let be h m the mth Hermite function. It is known that the ambiguity functions of Hermite 
functions are given by the Laguerre functions [34] (see for example also [5]). Obviously, the 
maximal eigen-value fulfills: 

A m ax(<5) > (h m ,Qh m ) = 7—7 / (h m ,S { ^h m )dfi = 7—7 / l m (n\n\ 2 )d[i (56) 

\ u I Ju \ u \ Ju 

where l m (t) =e~ t l 2 lJ${t) are the Laguerre functions and L,y are the 0th Laguerre polyno- 
mials. ■ 



E. Gaussian Signaling and the Corresponding Bounds 

The previous part of this section indicates that approximate eigen-functions have to be 
"Gaussian-like". Hence it makes sense to consider Gaussian signaling explicitely. For sim- 
plicity we do this for the time-frequency symmetric case g — 7 = 2^e~ 7Tt2 and n — 1. We 
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have the relation: 

(S^g)(x) = e -<2irX^e)+i(l-2c*)C(M+^), 7 ( x ) (57) 

if we let e := (l,i). According to (33) the error E p (fi) can be calculated as: 

EM = || [ dvV%>(v)e- i2 ^M . g . .)|| p (58) 

The function / : R 2 x R — > C is defined for a particular polarization a as: 

{ e -7r[2«7(^2 e )+i(l-2a)C(!^)+^ 2 ] _ e -7r[^-2iaC( l /,i/)] for C1 
(59) 
e -ir[2iTi(u,xe)+i(l-2a)C(v,u)+i%] _ j f Qr £2 

where we have used that the ambiguity function in polarization a is 

A£V) = e"W^ and 

Sa(^) : = v ■ v + 4ia£(z/, z/). The following Corollary contains the bounds specialized to the 
Gaussian case: 

Corollary 16 (Gaussian Bounds) For the case CI (k — 1) and for the case C2 (k = 2) in 
polarization a = it holds for any 1 < p < oo and 1 < q < oo that: 

EM < 32^ • ||A;(1- e-f») X X^ (60) 



I^Ti Ik 



where so(v) := (f, v). For g > 1 it follows from (60) a/50: 



^-<32^.k.(k\U\(l-(C,e-^)) 



1/ max((j',p) 



| V («) 

l^-H Iky 



(61) 



where C = Xu/\U\. 



Proof: We use the abbreviation Ai = e~ ns ° and A 2 = Re{e~^ SQ } as introduced in 
(27). Only for a = the case C2 provides an Euclidean distance measure in phase space. 
Equation (60) of the claim follows from Lemma 8 and from (39) of Lemma 9 together with 
Halloo < 2 1 1 1 1 oo = 32 1 / 4 . If q > 1 we can relate this further by (40) of Lemma 9 to weighted 
norms of ambiguity functions. Using the uniform bound C pq < k from Lemma 10 and the 
relation for Ri in (38) we get Gaussian integrals of the form (61) which can now be solved 



analytically for some cases. For example, if U is a centered disc of radius \/\U\/ir we get 
(C,e% s °) = l(2\U\/k) where l(x) = 2(1 - e~ x / 2 )/x. For a centered square of length \J\U\ 
we have instead l(x) = 2erf( ^Jnx/9>) 2 / x. ■ 
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V. Numerical Verification 

In this part we shall establish a spreading model with a finite number of random parameters. 
We shall need this model to verify numerically the bounds derived in this paper. Since several 
(iterated) integrals are involved which partially can only be computed numerically we have 
evaluate the achieved accuracy. We aim at computing E p (fi)/\\T,^\\ q up to a desired accuracy 
A. In our derivation we will assume that single definite integrals can be computed within a 
given predefined error (for example in using Simpson quadrature). 

A. Spreading Model with Finite Number of Parameters 

Let us consider a doubly-dispersive channel model with a finite number of fading pa- 
rameters Cfc, where k E 1? K and Z^- = {0...K — 1}. Each fading contribution has its 
own doubly-dispersive operation on the input signal, hence the model is different from the 
usual (distributional) models having a finite number of separated paths with fixed Doppler 
frequencies. The spreading function should be of the form: 

= c *Xu(^ - u(k + o))= ^ c kXi{v/u -k + o) (62) 
kei? K kez 2 K 

where Xu(y) = X[o,«](z/i)X[o,m](z/2) is the characteristic function of the square [0, it] x [0,it] =: 

[0,u] 2 and o = (|, |). Thus the latter is a disjoint partition of the square [0, Ku] 2 with 

area (Ku) 2 . In other words, if we fix the support of the spreading function to be \U\, then 

it follows for a K 2 -sampling of this area that u = ^J\U\/K. For such a model the q— 

norm of the spreading function as needed for the calculation of the ratio E p /\\t!^ || ? is: 

ll^w ll« = m2 ^' ? II c II'? where ||c|| 3 := (%2 k {ckl' 1 ) 1 ^ is simply the qth vector norm of the vector 

c = (. . . , Cfe, . . . ) G C R2 . Let us abbreviate I = l[k) = k + o. With (59) we get for the 

integrand in (58): 



/ ^\u)e-^^g(x)f(u,x)du= V c k ■ g(x) [ Xl (- - l)e~^^ f(u,x)du 



(63) 



F k (x) 

The approximate eigenstructure error reads now as E p (p) = lEfeez 2 c k m 9-Fk\\ P - For a = 1/2 
and case C2 the integral in F k (x) can be calculated explicitely. In general, however, F k [x) 
has to be computed numerically up to a certain accuracy 5 (it is a well-defined and definite 
integral). Thus, let the computed value F k (x) be such that pointwise \F k (x) — F k (x)\ < 5 
for all x and k. We would like to use F k (x) instead of F k (x) to compute the approximation 
E p (fi) on E p (fx). However we have to restrict the remaining indefinite integral over x to a 
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finite interval I := [-L, L\. With J we denote its complement in E, i.e. J := Observe 
from (59) that |/| < 2, hence |F fc | < 2u 2 in (63) and that for a Gaussian \\g\\P = 1/^p. If 
we choose txL > max( ^/log(2u 2 / 5) , 1) we have: 

2u 2 



\F k g ■ xj\\p < 2u 2 \\g ■ X j\\ P = 2n 2 erfc( v ^L) 1 / p < ~^e^ L 



Oil 2 T 9 TL>1 r0 

\9\\ P -^ry- p e- L2 < \\g\\ P -2u*.e- L2 <5\\g\ 



(64) 



For such a chosen L the integration with respect to x over the interval / = [-L, L] can be 
performed again within an accuracy of 5. This yields for the overall calculation error: 

\E p (ji) - em\ <s + Yl N " F k ] )9 ■ Xi\\ P + \\F k9 • Xj\\p] 



k 

< (l + 2||c||i- \\g\\ P )S= fl + 2||c||i -p~^) S 



(65) 



If we choose 5 = A ■ ■ (1 + 2||c||i -p 2 p) 1 (and L respectively) we can guarantee 

that the error on E p (jf) /\\T^\\ q is below A. 

Remark 2 (Interference Estimates for Statistical Models) Consider the following exam- 
ple: The transmitter sends the signal through the unknown channel 7i. Let us again for 
simplicity use the finite-parameter spreading model (62) for a support U of square shape. 
The receiver already knows the vector of fading parameters c for the spreading function 
of the channel, the pulse g and 7 and the time-frequency slot fi. The normalized q-norms 
c q = \\c\\ q ■ K~ 2 / q of the K 2 fading coefficients characterize the statistical model for the 
spreading such that ||S^|| ? = If/I 1 / 9 • c q . If the contribution of this particular slot fi is 
removed from the signal it remains e := TiS^j — X^S^g. Let us assume that the receiver 
expects another information in the span of the function f (for example f = S u g could be 
another slot v). The interference will be (f,e). Let be Af(p) = \\f\\ p > ■ 32 4 p . We have 

|{/,e)| < EM ■ \\f\\ p , < A f (p) ■ (\U\(1 - £ / 2 )) 1 /max( 9 ', P ) . m i/ q . Cq (66) 
With the assumption that \U\ < 1 we use |[/|V™*(?',p)+i/? < |[/| such that: 

\(f, e)\ < Af{p) • (1 - L 2 ) 1/max(c? ' p) • W\ ■ c q (67) 

This means, for different statistical models ( characterized by c q ) and functions f ( character- 
ized by || /||p' in the quantity Af(p)) we can characterize the amount of interference. 
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B. Numerical Experiments 

We will consider now the case where the coefficients Ck of the vector c £ C r2 are identical, 
independent and normal distributed which refers to the doubly-dispersive Rayleigh fading 
channel. The square shaped support U has a random size \U\ taken from a distribution 
uniformly on the interval [10~ 3 , 10~ 2 ] corresponding to values of the time-frequency spread 
relevant in mobile communication. Each realization of the fading factor c and u = \J\U\ 
parameterize via (62) a random spreading function £^ in a given polarization a which 
give itself rise to a random channel operator 7i by Lemma 1. On this random channel we 
have evaluated E p (ji) for Gaussian signaling as described previously in Section IV-E. For 
each realization we have taken /i uniformly from [— 5,5] 2 . We have calculated jV = 1000 
Monte Carlo (MC) runs for different values of p and q. For each run E p (p)/\\ || 9 has been 
computed (corresponding to one point in Fig.l and Fig. 2) up to an accuracy of A = 10~ 8 . 
The computed values E p (p) are compared to the uniform bound in (42) of Lemma 11 which 
depends only on the support and is valid for any normalized g and 7. Improved bounds are 
valid only for particular g and 7 like the Laguerre/Gauss (GL) bound from Theorem (7). Fig.l 
shows the case CI for p = q = 2, where we expect the most tight results. The GL bound 
improves the uniform estimates approximately by a factor of 10. However the computed MC 
values are still below this estimate by a factor of approximately two. The latter estimate 
degrades to a factor of approximately 10 for p = 3 and q = 3/2 as displayed in Fig. 2. 

VI. Conclusions 

In this paper we have considered doubly-dispersive channels with compactly supported 
spreading. We have shown to what level of approximation error a description as simple 
multiplication operators is valid. We have focused on two well known choices of such a 
description, i.e. the multiplication with the (generalized) Weyl symbol of the operator and 
the case of Wigner smoothing. We found that in both cases the approximation errors can 
be bounded by the size of the support of the spreading function. Our estimates improve 
and generalize recent results in this direction. Furthermore we have drawn the relation to 
localization operators and fidelity measures known from the theory of pulse shaping. Finally, 
we have verified our estimates using Monte Carlo methods with a precise control of the 
numerical uncertainties. 
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Fig. 1. Approximate Eigenstructure for the case CI, p = 2, q — 2: Verification of 1000 Monte Carlo runs with the uniform 
bound in Lemma 1 1 and the optimized Laguerre/Gauss bound of Theorem 7. 




Fig. 2. Approximate Eigenstructure for the case CI, p = 3, q — 1.5: Verification of 1000 Monte Carlo runs with the 
uniform bound in Lemma 11 and the optimized Laguerre/Gauss bound of Theorem 7. 
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Appendix 

A. Proof of Lemma 3 

The following proof is motivated by [6]. 
Proof: For each complex Hilbert space with ||o:||| = (x,x) the following inequality 

\\x-y\\l < ]|z||2- \\y\\l +2 \(y,x-y)[ (6g) 

(a) (b) 

holds. Now let x = HSf^ and y = L { ^{n)S^. Using (10) the following upper bounds 

(a) = ( 7) (s£»1CHSP - L% n + L% n - \L^) 7) 

+ |( 7 , (s^H*HS^ - L#n) 7)| 
+ | / ^ n (u)e- t2 ^^ ( A# (v) - l)dv\ 

+ ||E^ w n||i 

(b) = \L%(pL)\ ■ |(7, (SP*HSP - L#) 7)| 

JR 2 " 

will give the proposition. ■ 



< 


|r(«) 


I r(°)l 2 




|r(«) 
\ Aj -H*Ti 




< 


|r(«) 


1 r( a )| 2 



(69) 



B. Proof of Lemma 9 

Proof: Firstly - note that Holder's inequality for the index pair (l,oo) gives V£ < 

p-2 2 

V^~ 2 ■ V 2 2 with equality for p = 2; and in turn \\V p \\q> < \\Voo p ■V 2 p \\ q >. For q > 1 we can 
rewrite this and use again Holders inequality. We get: 

ii^ii* *= wv^ 21 ■ yf ii^' < mi? • \\y?\\% ( = imi? • <; P aw 

For q — 1 (i.e. </ = oo) we obtain rhs of the last equation directly: 

Halloo = \\vtf ■ v} iu < ||y||? • i?L /p (7i) 
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which proves (39) of this lemma. From the definition of R in (36) it is obvious that the 
minimum of the bounds is taken at B(U) = A$(U) which is provided by CI. Because 
there it holds always equality for p = 2 this is also the optimizer for || V-jH^ for any q. From 
(70) we get further for p < q' < oo: 

R% = Rli p \\R/Roo ■ xullX < RWR\ /q ' (72) 

because in this case (R^/R^y/v < R{y)/Roo for all v e U. For q' = p equality occurs 
in the last inequality. This proves (39) of this lemma for q' > p. For q' < p we use the 
concavity of R q '/ p , i.e. we proceed instead as follows: 

R%=(\U\-\\R q ' ,p Xu/\U\ ||i) W 

V J (73) 

< \U\W ■ \\R X u/\U\ \\Y P < \U\& ■ R\ /p 

The bounds (72) and (73) agree for q' = p and are tight for q = p = 2. ■ 
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